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Abstract. A study is presented of the central structure 
and kinematics of the galaxy NGC 6951, by means of 
broad band B' UK images and high resolution high dis- 
persion longslit spectroscopy, together with archival HST 
WFPC2 V and NICM0S2 J and H images. We find that 
there is little ongoing star formation inside the bar dom- 
inated region of the galaxy, except for the circumnuclear 
ring at 5 arcsec radius. There is some evidence that this 
star formation occurs in two modes, in bursts and con- 
tinuously, along the ring and inwards, towards the nu- 
cleus. The equivalent width of the Can triplet absorption 
lines show that, in the metal rich central region, the con- 
tinuum is dominated by a population of red supergiants, 
while red giants dominate outside. The gaseous kinematics 
along three slit position angles, and the comparison with 
the stellar kinematics, suggest the existence of a hierarchy 
of disks within disks, whose dynamics are decoupled at 
the two inner Linblad resonances (ILR), that we find to 
be located at 180 pc and at 1100 pc. This is supported 
by the structure seen in the high resolution HST images. 
The nucleus is spatially resolved in the emission line ra- 
tio [Nil] /Ha, and in the FWHM of the emission lines, 
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within a radius of 1.5", just inside the innermost ILR. 
Outside the ilLR, the stellar CaT velocity profile is re- 
solved into two different components, associated with the 
bar and the disk. Several results indicate that this is a 
dynamically old system: the little ongoing star formation 
inside the bar dominated part of the galaxy, the very large 
relative amount of molecular to total mass within the inner 
6 arcsec radius, ^ 25%, and the geometry of the circum- 
nuclear ring that leads the stellar bar at a position angle 
greater than 90°. It is thus possible that a nuclear bar 
has existed in NGC 6951 that drove the gas towards the 
nucleus, as in the bars within bars scenario, but that this 
bar has already dissolved by the gas accumulated within 
the circumnuclear region. We discuss the possibility that 
the kinematical component inside the ilLR could be due 
to a nuclear outflow produced by the combined effects of 
SN and SN remnants, or to a nuclear disk, as in the disk 
within disk scenario that we propose for the fueling of the 
AGN in NGC 6951. 

Key words: galaxies: active - galaxies: individual: NGC 
6951 - galaxies: kinematics and dynamics - galaxies: nu- 
clei - galaxies: spiral - galaxies: structure 



1. Introduction 

The centers of spiral galaxies are frequent sites of activ- 
ity, marked by the presence of intense star formation (SF) 
and/or an active galactic nucleus (AGN). This activity 
needs to be fueled with a supply of gas, whose reservoir 
can be provided by the disk of spiral galaxies. An efficient 
way for this gas to lose its angular momentum is pro- 
vided by the existence of non-axisymmetric components 
of the galactic potential, together with a strong central 
gravitational potential. The latter is necessary to explain 
that AGNs, irrespective of the presence of morphological 
perturbations, appear preferentially in early type spirals, 
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with the frequency peak at Sb. On the other hand, galax- 
ies dominated by disk star formation occur preferentially 
in later types with a distribution peaking at Scd (Moles 
et al. 1995; Ho et al. 1997a and references therein). 

The existence of gas to fuel the circumnuclear activ- 
ity is necessary but not sufficient (e.g., Moles et al. 1995; 
Maoiolino et al. 1997; Mulchaey & Regan 1997). The right 
dynamical and physical conditions must exist for this gas 
to be used effectively in either infalling to the nucleus 
proper, and feeding the AGN or nuclear starburst, or col- 
lapsing by self-gravity in the circumnuclear region in the 
form of intense SF. There are examples of galaxies where 
the gas is known to be available in the central regions, 
but where the onset of SF has yet to occur. For example, 
NGC 4151, which harbours a Seyfert type 1 nucleus, has 
a circumnuclear ring of material discovered in the form of 
dust extinction, and has a large column of neutral hydro- 
gen towards the nucleus, but there is no significant star 
formation occurring at present (Robinson et al. 1994; Vila- 
Vilaro et al. 1995; Mundell et al. 1995). However, in the 
circumnuclear ring of star formation in the target galaxy 
of this paper, NGC 6951, Kohno et al. (1999) find that al- 
though the dynamical shocks provide the mechanism for 
the accumulation of molecular gas along the ring, this may 
not be the mechanism responsible for the star formation, 
but rather the gravitational instability of the clouds thus 
formed is the driving force for the onset of star formation. 

To understand what are the conditions and mecha- 
nisms for the onset of nuclear activity in spiral galaxies, 
a detailed characterisation of morphological and kinemat- 
ical components in galaxies of different morphology and 
activity level is needed. The DEGAS collaboration was 
established to address such a detailed study in a small 
sample of spirals with an AGN, selected as isolated to 
avoid the possible external contribution to the assymme- 
trics in the potential. We plan to procede with a parallel 
analysis of a sample of isolated non-AGN spirals, to look 
for differences and/or similarities and their implications. 
IR imaging data of the complete sample are presented 
and analysed in Marquez et al. (1999a, 1999b). Other re- 
cent IR imaging studies of active and control samples in- 
clude McLeod & Rieke (1995), Mulchaey et al. (1997), and 
Peletier et al. (1998). 

In this paper we present the first results for the 
SAB(rs)bc galaxy NGC 6951. There are a number of stud- 
ies of this galaxy that include broad and narrow band opti- 
cal and infrared images (Buta & Crocker 1993; Marquez & 
Moles 1993, hereafter MM93; Earth et al. 1995; Wozniak 
et al. 1995; Elmegreen et al. 1996; Friedh et al. 1996; Rozas 
et al. 1996a, b; Gonzalez Delgado et al. 1997; Gonzalez Del- 
gado & Perez 1997; Mulchaey et al. 1997), spectroscopy of 
the nuclear region (e.g., Boer & Schulz 1993; Filippenko & 
Sargent 1985; Muiioz-Tuiion et al. 1989; MM93; Ho et al. 
1995, 1997b; Ho et al. 1997c), radio interferometric maps 
(Vila et al. 1990; Saikia et al. 1994), and high resolution 
spectroscopic molecular maps (Kenney et al. 1992; Kohno 



et al. 1999). To better characterize the detailed kinemat- 
ical components and their relation to the morphological 
structures, we have obtained high resolution, high disper- 
sion spectroscopic observations of the gaseous and of the 
stellar components, and direct imaging in optical and in- 
frared bands. 

This paper is structured as follows. In section 2 we 
present the imaging and spectroscopic data. Sections 3 
and 4 detail the results from the analysis of the images and 
from the spectroscopy respectively, while in section 5 we 
discuss these results into a common picture of NGC 6951. 
Section 6 gives the summary and our main conclusions. 

2. Observations and data reduction 

The journal of observations, given in Table |l|, summa- 
rizes the spectroscopic and imaging observations carried 
out with different telescopes and instrumental setups, and 
those retrieved from the HST archives. 

2.1. Spectroscopy 

We observed NGC 6951 on the night of 9/10 August 1996 
with the ISIS double spectrograph attached to Cassegrain 
focus of the 4.2m William Herschel Telescope, and the 
gratings B1200 in the blue arm and R1200, simultane- 
ously, in the red arm. The observations include three wave- 
length ranges: one range around H/3 and one around Ha 
(both through the blue arm of ISIS), and a near infrared 
wavelength range around the Ca ii triplet absorption lines 
through the red arm of ISIS. At each slit position angle 
we took two 1800 s integrations in the red arm, each si- 
multaneous with one 1800 s in the H/3 and one 1800 s in 
the Ha ranges. 

The gratings provide a linear dispersion of 0.39 
A/pixel, and the wavelength ranges covered with the 
TEKl and TEK2 CCD chips used are 4757-5141 A, 6493- 
6881 A in the blue arm, and 8506-8882 A in the red arm. 
We used a slit width of 1 arcsec, that projects onto 2.1 pix- 
els (0.81 A) at the detector. The spatial sampling along 
the slit is 0.36 arcsec/pixel. NGC 6951 was observed at 
three slit position angles: the major axis at 138°, the mi- 
nor axis at 48°, and an intermediate position at 84°. 

The observing was performed under photometric at- 
mospheric conditions. At the WHT the seeing is measured 
continuously with a guide star; it was slightly variable dur- 
ing the observations of NGC 6951, with a median value of 
1.2 ± 0.1 arcsec. The phase of the moon was grey. 

The reduction of the spectroscopic data followed the 
standard steps of bias subtraction, flatfield correction, 
wavelength calibration with a CuNe lamp observed be- 
fore and after the target, atmospheric extinction correc- 
tion, and flux calibration using the spectroscopic stan- 
dards G24-9, LDS749B, BD284211 and BD254644 ob- 
served through an 8 arcsec wide slit. Cosmic ray hits were 
removed individually from the 2D frames. Finally the sky 
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background was subtracted from each frame using the out- 
ermost spatial pixels free from galaxy contribution. The 
two exposures in the CaT range were added together after 
checking that the spatial and spectral alignment was good 
to a small fraction of a pixel. FIGARO and IRAF^ were 
used throughout the calibration and analysis phases. 

2.2. Optical imaging 

Optical images through two continuum filters, B and /, 
were taken with the Brocam II camera on the Cassegrain 
focus of the 2.56m Nordic Optical Telescope during the 
nights of 8 and 9 October 1996. B' is a special purpose 
made filter with a square spectral response centered at 
4630 A and with a FWIIM=290 A, so that it does not 
include any bright emission line. The equivalent width of 
any emission line within the / filter is negligible. Imaging 
through emission line free filters is important to ascertain 
the true circumnuclear stellar structure, that can other- 
wise be significantly distorted with standard broad band 
filters, where the emission lines can make an important 
flux contribution (Vila-Vilaro et al. 1995). The detector 
was a 1024x1024 Tektronix chip with a sampling of 0.176 
arcsec/pixel and a 3x3 arcmin field of view. Weather con- 
ditions were photometric and the average FWHM of the 
seeing profile is 0.8 arcsec. The journal of observations is 
included in Table |l|. 

The reduction and calibration of the images has been 
carried out in IRAF. Bias and flat field corrections were 
done in a standard way. The task COSMICRAYS was used 
for the overall cosmic ray removal. 

The final image of NGC 6951 for a given filter is ob- 
tained by adding all the images for that filter properly 
aligned, sky subtracted and calibrated. The alignment, 
using cross-correlation techniques involving only bright 
stars, is made to an accuracy of 0.1 pixel. The sky back- 
ground is estimated by averaging the median flux in a 
few 40 X 40 pixel patches located outside the galactic disk. 
The 3cr background hmit magnitudes are 20.7 and 21.0 
mag arcsec"^ in B' and /, respectively. 

2.3. Infrared imaging 

Infrared images through the fllter K' (Ac=2.10 /i, 
width=0.34 fi) were taken with the MAGIC camera at- 
tached to the Cassegrain focus of the 3.5m telescope at 
the German-Spanish Observatory in Calar Alto on 1996 
September 26. The detector was a Rockwell 256x256 pixel 

^ The authors acknowledge the data analysis software FI- 
GARO provided by the Starlink Project which is run by 
CCLRC on behalf of PPARC. IRAF is the Image Analysis and 
Reduction Facility made available to the astronomical commu- 
nity by the National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in 
Astronomy (AURA), Inc., under contract with the U.S. Na- 
tional Science Foundation. 



NICM0S3 array, and the camera was set up in the high 
resolution mode of the f/10 focus, giving a samphng of 
0.32 arcsec/pixel and a field of view of 82"x82". The ob- 
servation procedure consisted in 3x3 mosaic patterns of 
the object and nearby sky, alternating the exposures. The 
mosaic on the galaxy has a 90% overlap between the com- 
ponent tiles. During the observations, the seeing was 1.0 
arcsec and the conditions photometric. Each galaxy frame 
was subtracted from a median sky frame and divided by 
the flat field. The nine frames were recentered using two 
or three field point sources and the nucleus, and subse- 
quently they were median averaged. The flux calibration 
was performed with observations of the Kitt Peak faint 
standard stars numbers 30, 31 and 32. 

NGC 6951 was observed on 1996 September 30, as one 
of the targets of the BARS collaboration, through the fil- 
ters J and Kshort (Ks) with the ARcetry Near-Infrared 
CAmera, ARNICA, attached to the 2.56m Nordic Opti- 
cal Telescope in La Palma. ARNICA uses a NICM0S3 
array that gives a pixel scale of 0.51 arcsec/pixel and a 
field of view of 130" x 130". We integrated 40 seconds on 
the object (dithered in 11 positions) and the same time 
on the sky, for a total integration time on the object of 
440 s in J and 880 s in Ks (with a 10" shift between the 
different images to facilitate bad pixel removal) . Since the 
size of the galaxy is much larger than the field of view, we 
repeated the process four times, with the central region 
of NGC 6951 in common, placed in one of the corners. 
This allowed us to obtain a final mosaic image with a to- 
tal integration time of 1760 s in J (3520 s in Ks) in the 
central 2' and 440 s in J (880 s in Ks) in the rest. Flat-fleld 
frames were obtained from the median of the sky frames 
taken along the night. Sky-subtracted and flat-fielded im- 
ages were then aligned by using a number of foreground 
stars. We observed 4 standard stars for the flux calibra- 
tion, leading to a photometric accuracy of 10%. For all 
the reduction and calibration steps we used the package 
sqiid within IRAF. We measure a FWHM for the stars 
in the final mosaic of 1.3". The 3(T limiting magnitudes 
are 18.8 in J and 19.5 in Ks. 

2.4- HST images 

We have retrieved archival HST images taken with the 
WFPC2 through the V filter (547M), with a single expo- 
sure of 300 s, and with the NICM0S2 camera through the 
J (HOW) and H (160W) fiffers, with exposures of 128 s 
each. The spatial sampling of the central PC image in the 
WFPC2 camera is 0.046", and 0.075" for the NIC2 camera 
of NICMOS. The useful field of view is ~12"xl2", which 
contains the circumnuclear region. 

3. Imaging results 

Fig. |l] shows the B' and Ks images of the galaxy. The 
most conspicuous features in both images are the bar, the 
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bulge and the spiral arms. The bar is better traced in Ks, 
since the distortions produced by the dust lanes are min- 
imized. The circumnuclear ring and dust lanes are more 
clearly traced in B' , where they appear as straight fea- 
tures from the bar edges towards the center, curving ap- 
proximately within the bulge radius. All the ring features 
reported by Buta & Crocker (1993) are visible: (pseudo) 
outer, (pseudo) inner and circumnuclear rings, with sizes 
in agreement with those given by them, and also with 
optical and/or near- infrared images previously published 
(MM93, Friedh et al. 1996). We note the isophote twist- 
ing in the central parts, from the bar to the circumnuclear 
ring. 

To parametrize the results of the photometric anal- 
ysis, we have obtained the isophotal profiles by using 
the task ellipse within stsdas . ctnalysis . isophote in 
IRAF, which uses the ellipse fitting method. All the stars 
in the frame were previously masked. The results are given 
in Fig. ^, where we plot isophotal magnitudes, ellipticity 
and position angle (PA) as a function of the square root of 
the radial distance for B' , I, J and Ks bands. They agree 
with previous results by Friedli et al. (1996), who observed 
NGC 6951 with a similar configuration, although our data 
extends 30 arcsec further out. Since the area covered by 
our optical images is smaller than for MM93, and con- 
sidering that the faintest isophotes we reach in NIR are 
still contaminated by the spiral arms, we have decided to 
use the disk major axis position angle and inclination as 
determined by MM93, i.e., = 138° and i = 42°. These 
values are compatible with the outermost ellipses as traced 
by the sharp-divided images, and are in good agreement 
with the kinematically derived position angle and inclina- 
tion, as we will see below. 

To quantify the bar extension we have computed the 
amplitude and phase of the Fourier transform applied to 
the deprojected image, as explained in Marquez & Moles 
(1996). The mode m=2 has a constant phase out to 45", 
that we adopt as the bar extension. 

The isophotal profiles indicate the presence of the cir- 
cumnuclear ring and bar components, that appear as an 
excess around 5" for the former and from 20" to 50" for 
the second. These excesses are clearer once the photomet- 
ric decomposition is made. We have fitted an exponential 
disk for the outermost radii, subtracted it from the profile, 
fitted an t^^^ for the resulting bulge, subtracted it from 
the profile, and iterated until convergence is achieved. The 
results from the decomposition should be taken with cau- 
tion since the bar contamination and the dust obscuration 
are very important. The results for the infrared image are 
presented in Marquez et al. (1999). 

B'-I and J-Ks full size colour images are shown in 
Fig. ^. In B'-I the reddest features trace the dust lanes, 
most outstanding in the western side of the bar, while the 
bluest colours trace the regions of active star formation, 
that occur in the spiral arms outside a 30 arcsec radius. 
Within the part of the galaxy dominated by the bar struc- 



ture there is very little ongoing star formation, except in 
the circumnuclear ring that is delineated by a number of 
blue knots. A comparison with the Ha image (fig. 5 in 
MM93; fig. 10 in Gonzalez Delgado et al. 1997) shows 
that, for most regions, the Ha flux correlates with the 
colour B'-I, however the bluest large region at 46 arcsec 
along position angle 20° does not have a correspondingly 
bright Ha emission, indicating that it may be a star form- 
ing region in the post-nebular phase. A detailed analysis 
of the stellar populations in NGC 6951 will be presented 
in a future paper. In the two colour maps, bulge and bar 
colours are very similar. 

The highest contrast in the J-Ks image occurs around 
the western side of the circumnuclear ring, and along the 
western dust lane close to the ring. To look at these in 
more detail, Fig. ^ shows B'-K' and I-K' colour maps, 
where is the Calar Alto image, that has better sam- 
pling and seeing. Along the dust lanes, the colour con- 
trast is higher in B'-IC, that results from a larger dust 
absorption in B' than in /. For the reddest regions inside 
the circumnuclear ring, the situation is the opposite, with 
redder colours in I-IC , in particular, close to the two ends 
of the ring major axis (NW-W and SE-S), and the SW 
zone where the ring connects with the bar. This higher 
contrast in I-lC is similar to J-Ks in Fig. ^. The smaller 
contrast of the dust in B'-K can be explained by an ex- 
cess of produced by the ongoing star formation in the 
ring. Thus the colour maps can be qualitatively under- 
stood in terms of the combined effects of dust obscuration 
and stellar population. 

The HST WFPC2 V image has been published and 
analysed by Barth et al. (1995), that study the distribu- 
tion of super star clusters in the circumnuclear ring. The 
HST NICMOS images have not been published yet. In 
Fig. H we show the J-H colour image. The red features 
contouring the circumnuclear ring are visible in great de- 
tail, mainly delineated as two short spiral arms that trace 
the ends of the oval. Closer to the nucleus, it consists of 
a number of small knots and structures that delineate a 
multi-spiral like structure. 

4. Spectroscopic results 

The emission line spectrum 

The three slit orientations used are depicted in Fig. |[ 
that shows the Ha contours (from the data of Gonzalez 
Delgado et al. 1997) on a B'/K' grayscale colour map. 
The three position angles correspond to the major axis 
(PA=138°), the minor axis (PA=48°), and an interme- 
diate angle (PA=84°). The slits cross the circumnuclear 
ring, but only one of the bright Ha knots. 

We have measured the emission line intensities at ev- 
ery spatial increment along the slit, by means of fitting 
gaussian components with LONGSLIT (Wilkins & Axon 
1991). Fig. 1^ shows the spatial variations of the Ha line 
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intensity and the Ha/H/3 ratio for position angles (a) 48°, 
(b) 84° and (c) the major axis at 138°. The Balmer hnes 
are affected by underlying absorption components in the 
nucleus and in the ring Hli regions. A detailed analysis 
of the stellar population contributing to the Balmer ab- 
sorption is beyond the scope of this paper, but a compari- 
son of the spectra with evolutionary population synthesis 
profiles of the Balmer absorption line profiles (Gonzalez 
Delgado et al. 1999b), indicates that the Balmer emis- 
sion lines should be corrected by a line core equivalent 
width of 1 A in the nucleus and 0.5 A in the ring; this 
result is found through the comparison of the absorption 
wings detected in the Ha and H/3 lines and in the nearby 
metal absorptions with the predictions of the models. The 
Ha/H/3 ratio shown in Fig. ^ has been corrected in this 
manner. The Balmer emission line ratio maps the effect 
of extinction by the dust lanes; it is high in the nucleus 
(Ha/H/3 > 10 at the center), then becomes smaller in 
the transition region between the nucleus and the ring 
of Hii regions (Ha/H/3 ~ 4 at a distance of ±1.5 arc- 
sec), and increases again across the ring to reach values of 
Ha/up ~ 10 in the outer side of the ring. This redden- 
ing across the ring is clearly produced by the dust lanes, 
as seen in Fig. ^. The extinction that corresponds to the 
above values of the Balmer ratio is c(H/3)=1.6, 0.4, and 
1.6 mag for the nucleus, the transition region and the ring 
respectively. Thus, just outside the partially resolved nu- 
cleus, at ±1.5 arcsec, the extinction reaches a minimum 
and the Balmer ratio takes the case B value. The HST 
image in Fig. |^ shows indeed little dust in that interme- 
diate region. The Galactic extinction towards NGC 6951 
quoted in the NED IPAC database is 0.88 mag in B, cor- 
responding to E(_B-y)=0.22 or c(H/3)=0.31, i.e. most of 
the extinction that we measure in the intermediate re- 
gion between the nucleus and the circumnuclear ring is of 
Galactic origin and not intrinsic to NGC 6951. 

Fig. H shows the line ratio [N ii]6583/Hq; together with 
the Ha and [Nil] fluxes. For all three PA the ratio be- 
comes larger than unity within the inner ±2 arcsec, reach- 
ing values of 5 in the nucleus. This large ratio is due 
to a strong increase in the [Nil] flux and not to under- 
lying absorption in Ha. Indeed, the [Nil] takes normal 
H II region values in the ring, but increases sharply in- 
wards. This region of high [Nii]6583/Ha is significantly 
more extended than the nuclear point spread function, 
and thus it is spatially resolved. The [Sii]/Ha ratio (not 
shown) also shows a qualitatively similar behaviour, with 
an increase from [S ii]/Ha=0.18 in the ring Hii regions to 
[Sil]/Ha > 1.8 in the nucleus. The electron density, com- 
puted from the line ratio [S il]6717/6731, increases from 
low values, Ng ^ 300 cm~'^, in the ring Hii regions to 
values larger than 1000 cm~'^ within the inner region and 
the nucleus (cf. Table |). These [Nil] and [Sii] line ratios 
seem to indicate the presence of a shocked component, also 
supported by kinematical evidence, as explained below. 



Fig. H shows the H/3 and [O III] 5007 line fluxes along 
PA 138°. The nucleus has a high excitation. The [Olii] 
flux decreases sharply to reach a ratio [O iii]5007/H/3 < 1 
outside ±2 arcsec from the nucleus. Thus the ring Hii re- 
gions have very low excitation; in fact, we cannot measure 
[O III] 5007 pixel by pixel outside the nuclear component 
at any of the three position angles. In order to compute 
this excitation ratio for the ring Hii regions, we have ex- 
tracted three one-dimensional spectra at each position an- 
gle, integrating the nucleus and the two sides where the 
slit crosses the ring. Fig. |l^ shows these for PA=48°. The 
measured H/3 and [O iii] 5007 line fluxes and the ratio 
[O iii]5007/H/3 are given in Table |, as well as the ratio 
with the H/3 flux corrected for underlying absorption (an 
absorption equivalent width of 1 A in the nucleus and 0.5 
A in the Hii regions). 

4-2. The absorption calcium triplet spectrum 

We have measured the equivalent width of the two main 
CaT lines at AA8542, 8662A, according to the method de- 
scribed by Di'az, Terlevich & Terlevich (1989, hereafter 
DTT). The nuclear spectrum (central 0.72 arcsec) is plot- 
ted in Fig. ^ for reference; the three CaT lines, Mgl 
A8807, and two of the hydrogen Paschen lines in absorp- 
tion are identified. The weakness of Pal4 indicates a neg- 
ligible contamination of the CaT lines by Pal6, Pal5 and 
Pais, that fall in their red wings. The results of the pixel 
to pixel equivalent width measurements along the slit, for 
the three position angles, have been combined and are 
shown in a single Fig. |l^. The lefthand side of this fig- 
ure, (-10 arcsec to the nucleus) corresponds to the NE-SE 
quadrant on the eastern side of the galaxy. The right- 
hand side of the figure (from the nucleus to 10 arcsec) 
corresponds to the NW-SW quadrant on the western side. 
There is a clear systematic trend in this combined figure, 
that can also be appreciated (albeit with a worse signal- 
to-noise ratio) along the three individual PA. Outside the 
circumnuclear ring, ew(CaT) takes a low value of ^ 5.5 
A; at and within the ring, 5 arcsec east and 4 arcsec west, 
ew(CaT) jumps to a higher value of 7 A; and closer in, 
within the partially resolved nucleus, ew(CaT) rises to a 
value between 8 and 9 A. These are three distinct regimes 
clearly present in the circumnuclear region. 

The interpretation of ew(CaT) as a function of a stellar 
population depends on the metallicity (DTT; see the most 
recent work by Garcia- Vargas, MoUa & Bressan, 1998, 
and references therein). We can have an indication of the 
metallicity by measuring the ew(Mgl). According to DTT 
this line is sensitive to metallicity and to effective tem- 
perature, but not to gravity. In the nucleus of NGC 6951, 
ew(MgI)> 0.8 A, that indicates (see fig. 8 of DTT) the 
production by cool (Tcff < 4700 K), high metallicity stars 
(solar or higher). For this metallicity range, ew(CaT) is 
a function mainly of gravity and so we can now interpret 
the measurements in terms of stellar populations. 



6 



Enrique Perez et al.; Structure and kinematics of NGC 6951 



The lowest values of ew(CaT), present mainly outside 
the ring, are due to an old stellar population of giants, 
with an age around 0.5 to 1 Gyr. The highest nuclear val- 
ues, ew(CaT)> 8 A indicate the presence of an important 
population of red supergiant stars (RSG), that dominate 
the luminosity at these wavelengths, with an age around 
10 to 20 Myr. In the circumnuclear region, ew(CaT) takes 
intermediate values; it is significantly higher than in the 
outer parts, but the actual value belongs to a range that 
can be explained by either giants or RSG. Our interpre- 
tation is as follows. The measured equivalent widths in 
the circumnuclear and nuclear regions are a lower limit 
to the actual values, for two reasons. First, these equiv- 
alent widths are diluted by the continuum of the young 
stellar population responsible for the ionization of the gas; 
we cannot precisely gauge with our present data set how 
important is this dilution, for which we would need spec- 
troscopy in a longer wavelength range, including the ul- 
traviolet, but we do see the blue stellar knots so that some 
dilution by them must be taking effect. Second, at least 
outside a radius of 2 arcsec, we can distinguish two distinct 
stellar kinematical components; although we cannot mea- 
sure independently their contribution to the luminosity or 
to the ew(CaT) , it is at least possible that one of these two 
distinct populations (one that can be identified with the 
bar population) contributes to dilute the other population 
(whose kinematics is dominated by rotation, see below). 
Thus, if dilution is taken into account it is likely that the 
circumnuclear values of ew(CaT)'^ 7 A are in fact larger, 
so that they also reflect a population of RSG within the 
ring. 

4-3. Gas kinematics 

We have measured the systemic velocity in Ha by integrat- 
ing all the spatial increments along the slit corresponding 
to the disk emission (i.e. excluding the nucleus and the cir- 
cumnuclear ring of Hii regions), and fitting the two peaks 
of the resulting integrated profile. The mean velocity of 
these two peaks is 1417±4 km s~^; this corresponds to 
a scale of 92 pc/arcsec, for Ho=75 km s~^ Mpc~^. The 
velocity curves obtained from the pixel to pixel measure- 
ments of the Ha emission line are shown in Fig. |l^; the 
Ha line fiux is also plotted for reference. Allowing for the 
different spatial and spectral resolutions, these velocity 
curves agree very well with those presented by MM93. In 
the following we analyze the detailed and the global kine- 
matic structures found along the three position angles. 

First we confirm that PA=138° presents the largest 
velocity amplitude, followed by PA=84° and PA=48°, in- 
dicating that 138° is closest to the major axis. At first 
glance, there seems to be much velocity structure at 
PA=48° that does not correspond to a kinematic minor 
axis, where we would not expect a systematic residual ve- 
locity curve. On a more detailed analysis, this velocity 
structure present along PA=48° is related to local mor- 



phological structure, and we can see that the outer parts 
of the disk velocity correspond to the systemic velocity, as 
indicated by a horizontal dotted line, 1417 km s^^. Thus 
we adopt 48° and 138° as the kinematic minor and major 
axes respectively. 

For this value of the kinematic major axis and an in- 
clination angle of 42° (MM93), the total dynamical mass 
enclosed within the circumnuclear region (a radius of 6 
arcsec, equivalent to 644 pc), is M = rv^/G — 5.8 x 10^ 
M0. Kohno et al. (1999) report a molecular gas mass of 
1.4 X 10^ Mq within a 6 arcsec radius; this amounts to 
25% of the total dynamical mass. 

Much detailed structure is present in the velocity 
curves at the different PA. For example, at PA=48° we 
detect the 50 km streaming into the main NE arm, 
between -56.5 arcsec and -53 arcsec (labelled A in Fig. |l3|). 
Also at this PA we detect the 50 km s^^ streaming into 
the bar region at -36 to -20 arcsec in the NE (labelled 
B), and at 32 to 23 arcsec to the SW (labelled B'). At 
PA=84°, the 50 km s~^ streaming into the bar is seen at 
the eastern edge, between 40 and 30 arcsec (labelled C), 
and at the western edge of the bar between 39 and 32 
arcsec (labelled C'). 

When deprojected for an inclination of 42° and a kine- 
matic major axis of 138°, the velocity curves at PA 84° 
and at PA 138° generally match in the outer disk regions, 
except for local distorsions produced by the spiral arms. 
However, the circumnuclear region shows a different pic- 
ture. One important feature of the velocity curve at PA 
48° is the existence of an apparent counter-rotation of the 
gas within the circumnuclear region with respect to the gas 
in the disk at the same PA. Indeed, when we look into the 
±8 arcsec circumnuclear velocity gradients at the three po- 
sition angles more closely (insets in Fig. |l3|), we see three 
features. First, if we deproject the circumnuclear veloc- 
ity curves using the same two parameters as for the main 
disk, i — 42° and (f) = 138°, the deprojected amplitude at 
PA=84° becomes significantly larger than at PA=138°; 
this would imply that the circumnuclear region requires a 
different set of deprojection angles. Second, there is a con- 
tinuous change of slope between the three PA, fiattening 
from PA=138° to PA=84° to PA=48°. If the circumnu- 
clear velocity curve along PA=48° is also produced by 
disk rotation, then the kinematic axis of this circumnu- 
clear rotation must be different from the kinematic axis 
of the main disk. Under this assumption of disk rotation, 
we can compute a kinematic major axis for this circumnu- 
clear region. We obtain that in this case the major axis of 
the circumnuclear region would correspond to 118°; with 
this circumnuclear major axis and for the same|^ inclina- 
tion of 42°, the three deprojected circumnuclear velocity 

^ If this circumnuclear region is a decoupled gaseous disk, it 
could have a different inclination to the line of sight that the 
main galaxy disk, but we will assume here that this is not the 
case and that the inclination is the same. 
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curves agree quite well. Third, at PA=84°, the circumnu- 
clear velocity amplitude, 200 km s~^, is significantly larger 
than the galaxy main disk velocity amplitude at the same 
PA, 130 km s~^. These three facts, namely, the appar- 
ent nuclear counter-rotation at PA=48° and the greater 
velocity amplitude at PA=84° (both, with respect to the 
disk amplitude at PA=84°, and with respect to the de- 
projected circumnuclear amphtud at PA=138°), argue for 
a gas dynamics in the circumnuclear region of NGC 6951 
that is decoupled from the dynamics of the main body of 
the galaxy. We shall return to this point in the discussion 
section. 

The FWHM of Ha and [Nil] are plotted in Fig. |l|. 
The values shown have been corrected for an instrumen- 
tal resolution of 43.5 km s~^. The influence of the nuclear 
velocity gradient on the widths is negligible; the gradient 
is 44 km s~^ arcsec"'^, that amounts to a quadratic cor- 
rection on the nuclear widths of only 2 km s^^. The Ha 
line flux is shown (dotted) for reference. The velocity dis- 
persion in the circumnuclear region is very high, in the Ha 
line there is a central plateau of 180 km s~^ in the nuclear 
±1.2 arcsec, it then decreases down to a minimum of 70 
km that is reached at the maximum of Ha emission 
in the Hii regions. After this local minimum, the veloc- 
ity dispersion increases again to 120 km and above 
in the two quadrants NE-NW and SE-SW across the ring 
of Hii regions. This trend is consistent with the general 
picture of the dynamics deduced from our data and from 
the molecular content (Kohno et al. 1999), that will be 
developed in the discussion section. The behaviour of the 
velocity dispersion for the [N ii] A6583 line is qualitatively 
similar to but systematically higher than that of Ha by a 
factor 1.2 in the inner ±3 arcsec, reaching values of up to 
220 km s~^ in the nucleus. These emission line widths in- 
dicate dynamically hot nuclear and circumnuclear regions. 
The kinematic data and the results from the emission line 
ratios and density values presented in section 4.1 seem to 
imply the existence of shocked material. 

The width of Ha is not resolved further out in the 
disk of the galaxy, except across the spiral arms where it 
reaches values of 25 to 45 km s~^. 

4-. 4- Stellar kinematics 

We have measured the velocity curve from the two main 
lines of the Can absorption triplet at AA8542, 8662 A, by 
means of cross-correlating the galaxy frames with seven 
spectra of five different K giant stars, observed during 
the same night and with the same setup. These stars 
are HD132737, HD171232, HD208817, HD198858, and 
HD1918. The cross-correlation is performed with a 2D pro- 
cedure we have implemented that makes use of the IRAF 
function CROSSCORR. The position and amplitude of 
the peak of the cross-correlation function is measured, and 
the resulting velocity curves are weight averaged at each 
position angle. The signal-to-noise ratio of the data allows 



measurements only within a 20 arcsec radius at PA=84° 
and 48°, and 10 arcsec at PA=138°. The resulting curves 
are shown in Fig. |l5| (expanded for the central region in 
the insets) for the three position angles, together with the 
corresponding Ha emission velocity curves for comparison. 

First, we inspect the stellar velocity curves. Inwards of 
the circumnuclear ring, in the ±2 arcsec radius resolved 
nuclear zone, there seems to be a normal stellar rotation 
curve, with a measured gradient of 41 and 26 km s~^ 
arcsec- 1 at PA=138° and 84°, and flat at PA=48°. Out- 
side this radius it is possible to partially resolve the ab- 
sorption line profile into two components, that are seen 
in the cross-correlation function by corresponding local 
peaks; in most spatial increments one component fully 
dominates the profile, while in some increments the posi- 
tion of the two components can be more easily measured. 
The velocities corresponding to these two components are 
plotted with filled and open symbols in the figure. On ei- 
ther side of the nucleus, these two stellar components have 
their velocity of the same sign, although their amplitude 
difference reaches 50 km s^"'^. 

We now look into the comparison with the gas Ha ve- 
locity curves. Where only one stellar component is mea- 
sured, i.e. within the nuclear 2 arcsec radius, the velocity 
curves of the gas and of the stars at PA=84° and 138° 
follow each other approximately. This is not the case at 
PA=48°, where the stellar velocity curve is flat, while the 
Ha curve presents a slope (c.f. previous section 4.3) with 
a peculiar structure in the central ±1.5 arcsec. 

At PA=84° the nuclear ±1.5 arcsec part of the gaseous 
and stellar velocity curves both follow the same similar 
pattern; however, when the two stellar components can 
be kinematically resolved, neither of them two match the 
gas kinematics. One component is very fiat with a mean 
velocity of 1440 km from 4.5 arcsec eastwards of the 
nucleus, and a velocity of 1400 km s^^ from 5 arcsec west- 
wards of the nucleus. This stellar component matches the 
velocity of the gas entrained in the bar both at 30 arc- 
sec east (at 1440 km s^^) and at 30 arcsec west (at 1400 
km s^^; c.f. previous section 4.2 and Fig. p^ ). We shall 
refer to this component as the stellar bar component. The 
other stellar component is steeper but not as much as the 
gas component between ±3 arcsec and ±7 arcsec, where 
the gas presents a significantly larger rotation amplitude 
than the stellar disk rotation at the same projected dis- 
tance (a similar case to that reported for the barred galaxy 
NGC 6701 by Marquez et al. 1996, cf. their fig. 12). The 
extrapolation of this stellar component meets the gas ro- 
tation in the disk at 40 arcsec east (at 1490 km s""'^) and 
at 40 arcsec west (at 1330 km s""'^); we shall refer to this 
component as the stellar disk component. 

At PA=138° the stellar velocity presents a similar be- 
haviour. Only one component is apparent within a 3 arc- 
sec radius; this generally follows the gas velocity, although 
with less structure in the nucleus. Outwards of this radius 
two stellar velocity components can be distinguished. At 
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the south-east the steeper of the two components reaches 
a maximum velocity of 1553 km s~^, that extrapolates to 
match the disk gas velocity. 

The stars and the gas kinematics are most different at 
PA=48°. The gas presents an apparent counter-rotating 
component within a 6 arcsec radius with respect to its 
outer disk rotation, and an additional peculiar velocity 
structure in the nuclear ±1.5 arcsec. The stars show a 
slowly rising kinematics from north-east to south-west, 
compatible with the disk rotation, and there is no indi- 
cation of any peculiar velocity structure in the nucleus 
similar to the gaseous one. This peculiar structure in the 
gas velocity along the minor axis could be due to an off- 
centered slit, to a gaseous nuclear outflow, or to a structure 
related to the ilLR. 

5. Discussion 

We shall now discuss the results obtained in the previous 
sections and try to put them into a coherent picture of the 
morphology and kinematics of the circumnuclear region in 
NGC 6951. 

5.1. The level of activity in the nucleus: Seyfert 2 or 
LINER? 

The active nucleus in NGC 6951 has been classified in the 
literature as a Seyfert 2 (e.g., Boer & Schulz 1993; Ho et al. 
1995, 1997b) or as a LINER (Fihppenko & Sargent 1985; 
Mufioz-Tufion et al. 1989; MM93). Of the different classi- 
fications, that of Ho et al. (1997b) is the most accurately 
performed, because of the data quality and the correction 
for underlying absorption, that affects somewhat to Ha, 
and more to H/?. The classification is based on the three 
Veilleux & Osterbrock (1987) diagrams [O iii]5007/H/3 
vs. [Nii]6583/Ha, [O iii]5007/H/3 vs. [Sii]6716-h6731/Ha, 
and [Oiii]5007/H/3 vs. [Oi]6300/Ha. 

We have extracted the central 1.0 x 1.2 arcsec as 
the most representative spectrum of the nucleus, and we 
measure [O iii]5007/H/3=2.59 (5.89), [Nii]6583/Ha=7.00 
(3.72), and [S ii]6716-h6731/Ha=3.06 (1.63). For the val- 
ues in parentheses, H/3 and Ha have been corrected for 
an underlying 1 A of equivalent width in absorption (as 
discussed in section 4.1 above); these are the values used 
in the diagrams. Further, recently we have had access 
to a lower resolution spectrum, from where we can mea- 
sure [Oi]6300/Ha >0.35. Using these values in the dia- 
grams, the nucleus of NGC 6951 is located in the very 
high excitation end of the region of these diagrams filled 
by the LINER points. This explains the uncertain classifi- 
cation, even when the correction for underlying absorption 
is taken into account. In fact, because the most critical 
value is [Oiii]5007/H/3 which is most sensitive to this ab- 
sorption correction, it makes the precise definite classifica- 
tion somewhat arbitrary. Ho et al. (1997b) classification 
criteria (cf. their table 5) consider a horizontal bound- 



ary between Seyfert 2 and LINERs at [O iii]5007/H/? >3; 
this puts NGC 6951 in their Seyfert 2 zone. Other au- 
thors consider that the boundary between the two classes 
changes with the value of the abscissa (e.g. Gongalves et 
al. 1999), and this would put NGC 6951 in the higher ex- 
citation end of the LINER zone. In either case, the value 
of [Nii]6583/Ha=3.72 is also very high, that may suggest 
an overabundance of nitrogen; this will be discussed in 
a future work on the stellar populations and gas phase 
abundances. 

In summary, the active nucleus in NGC6951 can be 
considered as a transition object between a very high 
excitation LINER and a possible nitrogen overabundant 
Seyfert 2. 

5.2. Colour maps: dust and stellar populations 

According to numerical simulations and to observations 
of barred galaxies, the interaction of the bar and the disk 
creates mildly shocked zones that enhance the formation 
of Hii regions, mainly at the ends and along the leading 
edges of the bar. This also produces a net fiow of disk gas 
and dust towards the central regions of the galaxy, where 
the interaction with the bulge and nuclear dynamics of- 
ten enhances the circumnuclear star formation later in the 
evolution of these systems. In Fig. ^, the scarcity of star 
forming regions in the bar (both, at the ends and along 
the leading front of the bar), and the presence of the cir- 
cumnuclear ring of star formation, argues for a relatively 
old dynamical age of the bar and associated systems in 
NGC 6951. Other imaging and spectroscopic results give 
further support to this scenario. 

When we consider all the information available (cf. sec- 
tions 3 and 4), a pattern seems to emerge in the circum- 
nuclear region in which: 

(a) There is a widespread population of supergiant stars 
that dominate the light along and inwards of the circum- 
nuclear ring, between about 6 arcsec and 1.5 arcsec radius. 
These are most clearly seen in the K dominated J-K map, 
together with the large CaT equivalent widths measured 
in this same zone (cf. section 4.2). 

(b) The correspondence between the regions more lumi- 
nous in B' and in Ha is only good for the two brightest 
ones, but there are two regions bright in Ha with appar- 
ently no corresponding enhanced B' that must be under- 
stood in terms of extinction; and there is widespread en- 
hanced B' along and inwards of the ring with little or no 
corresponding Ha emission, that indicates star formation 
in the post-nebular phase. It is interesting to notice that 
5' is enhanced both locally in a knotty structure along 
the ring and more diffuse along and inwards of the ring; 
this might be interpreted as a signature of two modes of 
star formation coexisting here: continuous and bursting. It 
is possible that a continuous star formation process with 
ocasional bursts takes place in a region such as this one, 
where a more or less steady supply of gas has been going 
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on for a long time, as supported by the large fraction of 
molecular to total mass within the region (cf. section 4.3 
above) . 

(c) There is some indication of a time sequence in the star 
formation history along the ring, more apparent in the 
W-N quadrant. First, the strong K knot at position angle 
~250°. Second, two knots bright both in and _?sr with 
very little diffuse Ha emission at PA=270° and 300°; this 
is where Kohno et al. (1999) find the leading secondary 
maximum of HCN distribution, and also the location of 
the leading secondary maximum of the 6 cm radio emission 
(Saikia et al. 1994). Third, the Hii region to the north, 
PA=342°, that is brightest in B' , Ha, HCN and 6 cm. Fi- 
nally, the region just behind this one along the ring, partly 
associated with the CO peak from Kohno et al., and that 
might be predicted to be the next in the sequence of star 
formation. A similar, albeit somewhat less well defined, 
trend is observed in the south-eastern side. 

In conclusion, a qualitative look at the colours suggests 
that there may be two coexisting modes of star formation, 
continuous and bursting, and that the bursting mode may 
be sequential from the spearhead backwards along the en- 
trained material as it is shocked in the circumnuclear ring. 

Some of the ideas proposed here are only suggested 
by the wealth of information indicated by the data, and 
a detailed quantitative analysis of the stellar populations 
will be necessary to elucidate the validity of some of them. 
We are gathering spectroscopic data at other wavelengths 
to address this future analysis. 

5.3. Linear analysis of the kinematics 

By performing a standard linear analysis of the kinemat- 
ics (e.g. Binney & Tremaine 1987) it is possible to obtain 
the location of the main resonances, co-rotation (CR) and 
the inner Linblad resonances (ILRs). We use results from 
numerical simulations that place CR close to the bar semi- 
major axis (Athanassoula 1992), to infer a bar pattern 
speed of 3 km s~^ arcsec^^ (32 km s~^ kpc~^) ^. This 
implies the existence of two ILRs, at 2 arcsec (180 pc, in- 
ner ILR) and at 12 arcsec (1100 pc, outer ILR) ^. From 
OVRO CO observations of the central part of the galaxy, 

^ Computing the intersection between the bar pattern speed 
and the curve uj — Hi/2 (were u!{r) is the angular velocity and 
K(r) the epicyclic frequency) directly from the data points in 
the velocity curve is hampered by the dependence of n{r) on 
the derivative of the velocity, which implies that small errors 
and structure in the velocity as a function of radius translate 
into large fluctuations in K(r). In the case of NGC 6951 this 
problem is further exacerbated by the scarcity of data points 
in the region outside the circumnuclear ring, between 10 arcsec 
and 30 arcsec. 

The shallow radial dependence of ij-'(r) at this distance from 
the nucleus, and the steeper dependence closer to the nucleus, 
implies that a 10% or 20% uncertainty in the position of CR 
does not affect significantly the location of the ILRs. 



Kenney et al. (1992) obtain values of 180 pc for the ilLR 
and 460 pc for the oILR; this latter value is a lower limit 
derived from their upper limit to the bar pattern speed of 
66 km s^^ kpc^^, and thus is compatible with our results. 

5.4- Structure of the circumnuclear region 

The linear analysis of the kinematics indicates the exis- 
tence of a possible ilLR at 2 arcsec. Do we see a correlation 
between this dynamical feature and any morphological or 
spectroscopic properties? 

Morphologically, the Ha, colour, and HST optical im- 
ages all show structure within the inner 2 arcsec; i.e., it is 
not a point like nucleus. Friedli et al. (1996) have looked 
without success for the existence of a secondary bar, that 
could be possibly associated with the existence of the ilLR. 
There is no indication of such an inner bar either in our 
images. In fact, the HST optical image shows a clear in- 
ward spiraling structure down to 0.5 arcsec radius (46 pc). 

Spectroscopically, several diagnostics are spatially re- 
solved within the inner 2 arcsec radius, such as the 
[Nil] /Ha ratio (Fig. |), the FWHM of both Ha and [Nil] 
(Fig. [l4| ), and the equivalent with of the CaT lines in- 
dicating the presence of an important population of red 
supergiants in the nucleus (Fig. ^2|) . The gaseous and stel- 
lar velocity curves also show features within the ilLR: for 
the three PA, the two stellar components merge at the 
ilLR and are no longer distinguishable; the gas velocity 
curves at PA=84° and 138° flatten within a 1.5 arcsec ra- 
dius, while at PA=48° the gas velocity shows a peculiar 
structure within a 2 arcsec radius, that is not present in 
the stellar velocity. 

This peculiar structure along the minor axis could be 
due to an off-centered slit, to a gaseous nuclear outflow, 
or to a structure related to the ilLR. The first possibility 
is ruled out because this velocity structure is present only 
in the ionized gas but not in the stars, that implies that 
it is a physically existing feature, given that the two spec- 
tral regions were observed simultaneously and through the 
same slit. Of the two other possibilities, our data do not 
have sufficient spatial resolution to make a more definitive 
conclusion. However, we can make a few comments. A nu- 
clear outflow could be produced by the combined effect of 
the winds and supernova explosions of the recent past gen- 
eration of star formation indicated by the red supergiant 
stars we have detected. Saikia et al. (1994) give a total flux 
for the nuclear beam component of NGC 6951 at 6 cm of 
1.1 mjy; they argue that this non-thermal flux is most 
likely due to supernovae. This flux implies a luminosity of 
0.47x10^'^ watt Hz~^, that can be converted to a super- 
nova rate using the models of Colina & Perez-Olea (1992), 
to give 0.002 yr~^. A similar calculation with the 20 cm 
radio flux in the nuclear component given by Vila et al. 
(1990), that implies a luminosity of 1.03x10^° watt Hz"!, 
yields the same value for the SN rate. This is a low SN 
rate when compared with bright nuclear starbursts such as 
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NGC 7714, for which Gonzalez Delgado et al. (1999a) find 
0.07 yr""'^, but is qualitatively consistent with the faintness 
of the nucleus in Ha. If a nuclear outflow were the expla- 
nation for the kinematic features in the nucleus, we would 
only be seeing the approaching, blueshifted, side of the 
outflow with a speed of ~100 km s~^; the typical expan- 
sion velocity of a superwind bubble blown by a nuclear 
starburst in dwarf galaxies is of the order of 50 km 
(Marlowe et al. 1995); but is significantly larger, a few 
100 km s~^, for nuclear superwinds in Starburst galaxies 
(Gonzalez Delgado et al. 1998). The receding side would 
be presumably occulted by a combination of nuclear ob- 
scuration and spatial and spectral resolution. 

A third possibility suggested by the data is that the 
nuclear gas dynamics inside the ilLR is partially decoupled 
from that of the circumnuclear region which, in turn, is 
partially decoupled from that at larger scales. This would 
be a scenario of nested disks within disks, where a circum- 
nuclear disk accumulates mass from the outer main galaxy 
disk via the torques produced by the bar, and becomes suf- 
ficiently massive to decouple from the main disk dynamics 
at the oILR; part of this infall proceeds further into the 
nucleus, where a similar rotating structure (a torus or a 
nuclear disk such as those seen in HST images, Ferrarese 
& Ford 1999; Ford et al. 1998) decouples from the circum- 
nuclear disk at the ilLR. We notice that the orientation 
of the ring in NGC 6951 is leading with respect to the 
inflow of material from the bar, i.e. it is advanced from 
the perpendicular to the bar, and this relative orientation 
is seen in numerical simulations only when the system is 
relatively old (see, for example, figure 3 of Byrd et al. 
1994). In this context, it would be possible to interpret 
the change of the velocity curve in the three PA, within 
the inner 2 arcsec radius, as produced by rotation in a disk 
with a minor axis different from the 48° of the main disk, 
and different from the 30° of the circumnuclear disk (cf. 
section 4.3); it would have a value somewhere between 84° 
and 48°, so that at PA=48° we could observe the change 
of velocity slope with respect to that at 84°. 

The appearance of the HST images would support this 
scenario. Figure shows the HST Fand if sharp divided^ 
images; the 77 image shows a very uniform stellar light dis- 
tribution in the transition zone between the ring and the 
nucleus. However, the F image shows a structure spiralling 
inwards, delineated by the dust that can be traced right 
into the nuclear 0.5 arcsec. Thus, these high resolution 
images show both the lack of any inner bar-like struc- 
ture and the existing of a continuing spiralling into the 
nucleus, similar to that found in numerical simulations by 
Finer et al. (1995, cf. their fig. 4). A similar conclusion has 
been recently reached by Regan & Mulchaey (1999), who 
have looked for evidences of strong nuclear bars in WFPC2 



^ This process of enhancement consists in dividing the origi- 
nal image by a median filtered version, so as to enhance sharp 
features. 



and NICM0S2 images of a sample of Seyfert galaxies, and 
find evidence of the existence of this nuclear bar in only 
3 out of 12 galaxies studied, while the majority of the 
galaxies show a spiral morphology similar to what we find 
in NGC 6951. This does not necessarily implies that the 
'bars within bars' mechanism proposed to fuel the nucleus 
has not been at work in these galaxies (Shlosman et al. 
1989; Friedh & Martinet 1993); at least in NGC 6951 we 
have found several indications that point to a system that 
is dynamically old, in particular, the large percentage of 
molecular mass accumulated in the nucleus, ~ 25%, could 
have already occasioned the dissolution of a nuclear bar, 
and be working towards the desassembling of the large 
scale bar. 

This idea is inferred from extrapolations of three slit 
PA. While we understand that a complete 2D spectro- 
scopic mapping of this and other galaxies is required to 
confirm or otherwise dismiss it, it is also true that such a 
2D mapping at a spectral and spatial resolution equal to or 
better than our data is very difficult to obtain with present 
day instrumentation. Recent developments like SAURON 
(Miller et al. 1999) and similar instruments will have an 
important impact on such studies in the next few years. 

6. Summary and conclusions 

We have obtained broad band B' UK^ images and high res- 
olution high dispersion longslit spectroscopy for the ion- 
ized gas (around H/3 to [O ill] and Ha to [S ii]), and for the 
stellar populations (in the Can triplet lines around 8500 
A) of the galaxy NGC 6951. Together with archival HST 
VJH images, we analyse these data to study the central 
structure and kinematics and find that: 

— There is little star formation ongoing inside the bar 
dominated region of the galaxy, except for the circum- 
nuclear ring at 5 arcsec radius. There is some evidence 
that this star formation occurs in two modes, in bursts 
and continuously, along the ring and inwards towards 
the nucleus. 

— The nuclear spectrum shows both very high excita- 
tion and very strong [N ll] and [S ll] lines, making the 
classification of the AGN somewhat uncertain between 
a high excitation LINER and a possibly high nitrogen 
abundant Seyfert 2, depending mainly on the uncertain 
correction for underlying absorption in H/3. The elec- 
tron density varies between 300 cm^'^ and 1000 cm~^. 

— The equivalent width of the Can triplet absorption 
lines show that in the metal rich central region of this 
galaxy, within 5" radius, the continuum light is dom- 
inated by a population of red supergiant stars, while 
outside the circumnuclear ring the stellar population 
is that of giants. 

— We suggest that the gaseous and stellar kinematics 
along the three slit position angles can be interpreted 
as the existence of a hierachy of disks within disks, with 
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dynamics decoupled at the two ILRs, that we find to 
be located at 180 pc and at 1100 pc. This would be 
supported by the structure seen in the high resolution 
HST images. 

— The nucleus is partly resolved within a radius of 1.5" 
(just inside the ilLR) both in the emission line ratio 
[Nii]/Ha, and in the FWHM of the emission fines. 

— Outside the ilLR the stellar CaT velocity profile can 
be partly resolved into two different components that 
seem to be associated to the bar and to a disk. 

— We discuss the possibility that the kinematic compo- 
nent inside the ilLR could be due to a nuclear disk, as 
in the disk within disk scenario suggested above, or to 
a nuclear outflow produced by the combined effects of 
SN and SN remnants. 

— Several clues indicate that this is a dynamically old 
system: (i) there is little star formation ongoing inside 
the bar dominated part of the galaxy (except for the 
circumnuclear ring), (ii) the relative amount of molec- 
ular to total mass within the inner 6 arcsec radius is 
very large ~ 25%, and (iii) the geometry of the circum- 
nuclear ring leading at a position angle greater than 
90° from the stellar bar. It is thus possible that a nu- 
clear bar has existed in NGC 6951 that drove the gas 
towards the nucleus, as in the bars within bars sce- 
nario, but that this bar has already dissolved by the 
gas accumulated within the circumnuclear region. 
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Fig. 1. Full frame (a) and Ks (b) broad band images 
of NGC 6951. The images show clearly two spiral arms, 
the bar, dust lanes along the bar and the circumnuclear 
ring at 5 arcsec from the nucleus. The orientation in all 
the images is North up and East to the left. The scale is 
92 pc arcsec^^. The contour levels are (a) in B' 24.0 22.5 
21.5 20.5 19.5 19.0, (b) in Ks 21.0 20.0 19.5 19.0 18.5 18.0 
17.5 17.0 16.5 16 15.5 15.0. 

Fig. 2. Broad band photometric parameters as a func- 
tion of the root square of the distance to the nucleus. Sur- 
face brigthness (/x): upper pannels; ellipticity (e): middle 
pannels; and position angle (PA): bottom pannels. The 
parameters arc derived by fitting isophotal profiles to the 
images: ff (left pannels); / (middle left pannels); J (middle 
right pannels) and Ks (right pannels). 

Fig. 3. Full frame colour images, (a) Bf -I ami (b) J-Ks. 
B'-I traces very well the spiral arms and the circumnu- 
clear ring, where recent star formation is taking place. 
Dust lanes along the bar are also seen as dark features. 
Grayscales are displayed between B' -1=3.1 and 1.9, and 
J-Ks=-1.0 and -2.0. 

Fig. 4. colour images of the bar and inner regions of 
NGC 6951: (a) B'-K' and (b) I-K. The inner part of the 
ring is dominated by blue colours, except for the nucleus 
that appears as a red central knot. Grayscales are dis- 
played between B'-K =0.06 and 0.16, and 7-ii:'=0.55 and 
1.30. 

Fig. 5. J-H colour image of the central 12"xl2". These 
images were observed by HST+NICMOS. The spatial 
scale is 0.075 arcsec pixel" ^, and the grayscale has been 
displayed between J-i/=0.55 and 0.31 mag. The circum- 
nuclear ring is mainly delinated as two short spiral arms. 
In the inner part of the ring, the dust delineates a multi- 
spiral structure that ends in the nucleus. 

Fig. 6. B'-K colour image of the central 20x20 arcsec. 
Ha contours (from the data of Gonzalez Delgado et al. 
1997) are also plotted. Note that there is not a perfect 
correspondence between the blue knots in the B'-K' image 
and the nebular emission. The three slit positions for the 
spectroscopy data are also plotted. Grayscale is displayed 
between B'-K=Om and 0.16. 

Fig. 7. Ha/H/3 emission line ratios (dots) measured along 
the three slit PA. This ratio has been corrected for un- 
derlying stellar absorption as explained in the text. The 
theoretical ratio predicted by the case B recombination is 
plotted as a horizontal dashed line. The Ha flux along the 
slit is also shown as a full line. 
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Fig. 8. Ha (open circles + dashed line) and [N ii] (squares 
+ dashed) fluxes and their ratio [Nil] /Ha (filled circles + 
full line) measured along the three slit PA. Note that the 
large [Nil] /Ha ratio in the inner ±2 arcsec is due to a 
strong increase in the [N ll] flux. 

Fig. 9. H/3 (filled circles + dashed line) and [O ill] (open 
points + full line) fltixes along PA=138°. The [Olll] emis- 
sion is mainly concentrated in the inner ±2 arcsec, where 
the H/3 emission is very weak. 

Fig. 10. One dimensional spectrum, from 4750 A to 5150 
A, of the nucleus and of the two regions on either side of 
the nucleus along PA=48°. Note that these regions have 
low excitation, consistent with the gas being photoionized 
by stars. 

Fig. 11. Nuclear spectrum (the central 0.72 arcsec) in the 
Call triplet region. The most relevant lines are labelled. 
Mgl A8807 is very strong, indicating that the metallicity 
is solar. Paschen 14 is very weak, thus Pal3 and Pal5 
must contribute very little to the Call A8542 and A8662, 
respectively. 

Fig. 12. Equivalent width of the Can A8542 + A8662 

lines, measured using the line and continuum windows de- 
fined by Diaz et al. (1989). The equivalent width is > 7 
A in the nucleus and out to the circumnuclear ring, indi- 
cating the presence of red supergiant stars in the central 
5x5 arcsec of the galaxy. 

Fig. 13. Velocity curve (dots) obtained from the pixel to 
pixel measurements of the Ha flux (full line), derived fit- 
ting a gaussian to the profile of the emission line. Curves 
are plotted for the three slit positions. At PA=138° (pho- 
tometric major axis of the galaxy) the curve shows the 
maximiim velocity amplitud. At PA=48° the curve shows 
more structure than expected for the minor axis of the 
galaxy. The inset within each figure shows the circumnu- 
clear ±8 arcsec expanded; the three insets have the same 
scales in both distance and velocity axes, respectively. 

Fig. 14. FWHM of the emission lines Ha (filled circles) 
and [Nil] (empty squares). The velocity dispersion of the 
gas is maximum in the inner ±1.2 arcsec, and minimum 
at the position of the H II regions in the ring. The Ha flux 
is plotted as a dashed line. 



Fig. 15. Velocity curves from the two strongest CaT 

absorption lines (AA8542, 8662), measured by cross- 
correlation of the galaxy frames with spectra of K giant 
stars. The velocity curve of the ionized gas (dots) is shown 
for comparison with the stellar velocity curve (circles) . In 
the central region the CaT absorption lines arc resolved in 
two components, plotted as fllled and open circles. The in- 
sets show the central ±10 arcsec; the three have the same 
radial and velocity scales respectively, and the gas veloc- 
ity has been plotted as a full line. Notice how in the inner 
±2 arcsec, the stars and the ionized gas follow the same 
velocity pattern at PA=84° and 138°, but not along the 
minor axis, where the Ha curve presents more structure. 

Fig. 16. HST WFPC2 V (a) and NICMOS H (b) im- 
ages of the circumnuclear region. To show more clearly 
the morphological structure of the cicumnuclear region, 
the original images have been sharp enhanced by dividing 
them by the corresponding median filtered image. Note 
the spiralling structure delineated by dust in the V image 
down into the nucleus. 
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Table 1. Journal of observations 



telescope - PA (°) 


sampling^ 


wavelength (A) 




airmass 


WHT - 138 


0.36" 


4757-5141 


1800 


1.32 


WHT - 138 


0.36" 


6493-6881 
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1.36 


WHT - 138 


0.36" 


8506-8882 


3600 


1.34 


WHT - 84 


0.36" 


4757-5141 
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1.27 


WHT - 84 


0.36" 
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1.26 


WHT - 84 


0.36" 


8506-8882 


3600 


1.26 


WHT - 38 


0.36" 


4757-5141 
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WHT - 38 


0.36" 


6493-6881 
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1.26 


WHT - 38 


0.36" 


8506-8882 


3600 


1.26 


NOT 


0.18" 


& 


3000 


1.36 


NOT 


0.18" 


I 


660 


1.51 


NOT+ARNICA 


0.51" 


J 


1760 (440) 


1.25 


NOT+ARNICA 


0.51" 


Ks 


3520 (880) 


1.5 


3.5m CAHA 


0.32" 


IC 


300 


1.22 


HST - WFPC2 


0.046" 


V (547M) 


300 




HST - NICM0S2 


0.075" 


J (HOW) 


128 




HST - NICM0S2 


0.075" 


H (160W) 


128 





PA 


distance 


density 




(arcsec) 


cm~^ 


138° SE 


-4.5 


530 


138° SE 


-1.1 


1440 


138° nuc 


0.0 


1340 


138° NW 


1.0 


990 


138 


i..") 


290 



Spatial sampling in arcsec per pixel. The linear 



dispersion is 0.39 A pixel ^ for all the spectroscopic observations. 
Table 2. Electron density along the slit 



Table 3. Hisli excitation ratio [Oiii]/H/'? 



PA 


orientation 


extraction 


_ F(H/3) 


F(5007)/F(H/3) 


F(5007)/F(H/3) 






arcsec 


lO^'^" erg cm~^ 




corrected" 


138° 


SE 


-5.69,-2.82 


1.13 


0.08 


0.07 


138° 


nuc 


-1.40, 1.11 


0.20 


16.26 


5.51 
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NW 


2.54, 5.40 
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0.03 


0.03 
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-4.00,-1.80 


0.38 
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0.29 
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0.12 


48° 


NE 
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0.30 


0.25 


48° 


nuc 


-0.89, 1.25 


0.37 


16.57 


7.96 


48° 


SW 


1.97, 4.47 


0.23 


0.53 


0.39 



" H/3 corrected for an absorption equivalent width of 1 A in the nucleus and 0.5 A in the Hii regions. 
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